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ruling out the above models. Other than the Yamada 
plot, we note that some new experimental developments 
have also been reported. Among these are the detailed 
evolutions of the resonance peak and incommensurability 
with doping p] , p^] , and the dispersion connecting these 
two structures which has a downward curvature opposite 
in sign to a conventional magnon dispersion [ fl3|Jl4| ] . 

In view of these new experimental observations, we 
present in this paper a detailed investigation of the dop- 
ing and frequency dependences of the incommensura- 
bility and the resonance peak energy, and their rela- 
tionship based on the nested Fermi surface scenario. 
We start from the slave-boson approach to the two- 
dimensional (2D) t—t'—J model. The aim of choosing this 
approach is that we can incoprate the evolutions of the 
FS and the superconducting(SC) gap with dopings self 
consistently at the mean-field level and determines their 
values according to the well accepted input parameters 
t, t' , J. This enables us to carry out a quantitative study 
of the doping and frequency evolutions of spin fluctua- 
tions. In most previous studies |p|, ^0|Jl5| , a tight-binding 
dispersion and the magnitude of the SC gap which are 
inferred from related experimental data for a fixed dop- 
ing concentration are taken. Consequently, these inves- 
tigations are limited to a fixed hole doping. Our main 
results are: 1) We find that the equality between the 
low-energy incommensurability and doping density in the 
underdoped regime is reproduced. This may provide an 
alternative to its exclusive stripe-phase explanation. 2) 
The downward dispersion of the collective mode is re- 
produced and understood as arising from the Fermi sur- 
face geometry, and its doping dependence is presented 
to be compared with further experiments. 3) Below a 
crossover frequency, the incommensurate peak shifts from 
the (-7T, 7r + 5tt) direction to the diagonal direction due 
to the node-to-node excitations. Above the resonance 
frequency, the incommensurate peak occurs again. It is 
along the (ir, n + Sir) direction in the underdoped regime 
[x < 0.13) and along the diagonal direction in the opti- 
mally doped and overdoped regimes x > 0.14. 4) The res- 
onance energy has a linear relation with doping in the un- 
derdoped regime and saturates near optimal doping, then 
slightly decreases with doping in the overdoped regime. 
As a result, a proportionality between the low-energy in- 
commensurability and the resonance energy is obtained. 
Most of these results are consistent with experimental 



The doping and frequency evolutions of the incommensu- 
rate spin response and the resonance mode are studied based 
on the scenario of the Fermi surface topology. We use the 
slave-boson mean-field approach to the t — t' — J model and 
including the antiferromagnetic fluctuation correction in the 
random-phase approximation. We find that the equality be- 
tween the incommensurability and the hole concentration is 
reproduced at low frequencies in the underdoped regime. This 
equality observed in experiments was explained only based on 
the stripe model before. We also obtain the downward disper- 
sion for the spin response and predict its doping dependence 
for further experimental testing, as well as a proportional- 
ity between the low-energy incommensurability and the res- 
onance energy. Our results suggest a common origin for the 
incommensuration and the resonance peak based on the Fermi 
surface topology and the d-wave symmetry. 

PACS number: 74.25.Ha,75.40.Gb,74.72.-h 



The inelastic neutron scattering experiment plays an 
important role in the studies of the spin dynamics in 
high-T c superconductors. It can provide the momentum 
and frequency dependences of the dynamical spin sus- 
ceptibility. Over the past decade, the striking feature of 
the spin susceptibility observed in the momentum space 
is the incommensurate peak along the (tt, n + Sit) direc- 
tion and that in the frequency dependence is the 
resonance peak at the antiferromagnetic (AF) wavevector 
Q = (7r,7r) ^j^. These stimulate intensive experimental 
and theoretical studies. Recently, it is suggested that the 
presence of dynamic stripes is the origin of the observed 
incommensurate peak, at least for La2- £! ;Sr ;r Cu04 Q. 
One of its stronge supports comes from the experiment by 
Yamada Q et al. which shows an equality between the 
incommensurability and the doping concentration (Ya- 
mada plot) in the underdoped La2- a; Sr a ;Cu04. This 
equality follows naturally from the static stripe model, 
but has not been explained in any other way up to 
now j7j. However, it is now unclear how to explain the 
resonance peak based on this scenario. On the other 
hand, the gross features of the incommensurate spin fluc- 
tuations can also be explained based on the scenario 
of the nested Fermi surface(FS) |||| and the resonance 
peak is thought to be a collective spin excitation mode 
in this framework jl0]j|]. Obviously, more fine and de- 
tailed experimental data are helpful for selecting out or 
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data Jljlul and the others wait for further experimental 
testing. Therefore our result presents an alternative ex- 
planation of the incommensuration and resonance peak 
in the spin response based on the d-wave superconduc- 
tor with tight-binding dispersion. However, the abrupt 
saturation of the incommensurability near the optimal 
doping and in the overdoped regime which is observed 
both in I^-zSr^CuC^ and YBa2Cu306+ K is not 
found. Possible reason for this is discussed. 

In the slave-boson approach to the t — t' — J 
model P,p7|, the physical electron operators Cj CT are 
expressed by slave bosons bi carrying the charge and 
fermions fi a representing the spin; Ci a = bffi a . In 
the SC state, we consider the order parameters Ay =< 
fitfjl ~ filfjl >= ±Ao with the (i-wave symmetry and 
Xij = J2a < ftrfja >== Xo, in which bosons condense 
bi — »< bi >= t/x (x is the hole concentration). Then, 
the slave-boson mean- field Hamiltonian of the 2D t—t'—J 
model in the SC state is, 

ka k 

+2NJ'( X l + &l) 1 



where = —2{xt + J'xo)[cos(fc x ) + cos(k y )] — 

Axt' cos(k x ) cos(fcy) — /i is the dispersion for fermions, and 
A fe = 2J' A [cos(fe x ) - cos(fc a )], with J' = 3J/8. The 
normal state Fermi surface at doping x — 0.12 for this 
dispersion is shown in Fig.l. The mean-field parameters 
Xo, Ao and the chemical potential fi for different dopings 
are obtained from a self-consistent calculation fl7| ]. 

It has been shown j||l^,|ll| that the inclusion of the AF 
fluctuation correction is necessary to account for some of 
the spin and charge dynamics. Formally, this can be done 
by perturbing around the mean-field Hamiltonian, i.e., 
we write H = H m + H 1 with H the usual 2D t - 1' - J 
Hamiltonian, and treat H' as a perturbation. In prin- 
ciple, all fluctuations are included. However, different 
selection of subset diagrams may result in different kinds 
of fluctuations. For the spin fluctuation, the usual ran- 
dom phase approximation (RP A) selects a series of ring 
diagrams as shown in Ref. [p|, [l^|Jl9| ] . Then, the renormal- 
ized spin susceptibility is given by, 



x(q,w) 



xo(q,w) 



i + ^7gXo(q,^) 



(2) 



where 7 9 = cos(q x ) +cos(q y ) . xo is the bare spin suscepti- 
bility which comes from the fcrmionic bubble and is given 
by the usual BCS form @. As done before |,[|Jl|], we 
choose ij — 0.34 instead of 77 = 1 in order to set the AF in- 
stability at x = 0.02 which is the experimental observed 
value. 

Numerical calculations are performed by dividing the 
Brillouin zone into 1024x1024 lattices, with t = 2J,t' = 
— 0.45i and J = 0.13eV. The damping rate of the 
fermionic quasiparticles is simulated to be T = 0.004J, 



however we note that the incommensurability and the 
resonance energy is not subjected to the change of T when 
it is below 0.1 J. Therefore, the only adjusted parameter 
throughout the paper is ry which has been fixed above. 

In Fig. 2, the peak position of the imaginary part of the 
renormailized spin susceptibility is shown as a function of 
frequency, for doping concentrations x — 0.08,0.14 and 
0.20 which may correspond to the underdoped, optimally 
doped and overdoped regimes. At their resonance ener- 
gies (0.375 J, 0.53 J and 0.545 J, respectively), the q peaks 
are at S = Q, i.e. at (ir,ir) point. When the frequency 
is reduced from the resonance energy, the incommensu- 
rate peak occurs and is found to be in the (7r, tt + 5n) 
direction. In particular, a downward curvature for the 
peak dispersion is observed for various dopings. This 
is in agreement with experiments fl5| , |l4| . To understand 
its origin, we show the frequency dependences of the bare 
spin susceptiblity Im^o at 1 = 0.14 for different wavevec- 
tors q = (tt, 0.727r), (tt, 0.8it) and (7r,0.97r) in Figs. 3(a) 
and 3(b), and at q = (tt, 0.847r) for different dopings 
x = 0.08,0.14 and 0.20 in Figs.3(c) and 3(d). Firstly, 
it is important to point out that the bare spin suscepti- 
blity Im^o is incommensurate for all energies and dop- 
ings in the ranges considered here. Moreover, its peak 
is along the (n + Sn, tt) or (n, n + Sty) directions. It 
comes from the FS nesting effect as has been shown by 
Brinckmann and Lee In the current framework, the 
origin of the resonance peak is ascribed to a collective 
spin excitation corresponding to the pole of the renor- 
malized spin susceptibility, i.e. 1 + 77 J7qRcxq(Q, u>) = 0, 
and negligibly small Imxo(Q,^) [|8], |l0| , |l8| . At the com- 
mensurate wavevector Q = (7r,7r), there is a step-like 
rise in Im^o (shown as the soild line in Fig. 3(b)) at the 
threshold energy for particle- hole (p-h) excitations. Due 
to this step-like rise, a logarithmic singularity in Rexo 
occurs via the Kramers-Kroenig relation. This singular- 
ity shifts downward the collective mode energy and leads 
it to situate in the spin gap, so no damping is expected 
for the mode. For the energy band efc and the SC gap 
considered here, the initial and final states of the p-h 
excitations at the threshold energy (minmum excitation 
energy) lie on the Fermi surface. So, the threshold exci- 
tation with transition wavevector Q corresponds to the 
excitation A — > B as schmatically shown in Fig.l. The 
step-like rise is caused mainly by the van Hove(vH) singu- 
larity around (tt, 0) and the cZ-wave symmetry of the SC 
gap jl^jlj,|l^]. In this case, the collective spin excitation 
at (tt, 7r) which is determined by the pole condition in 
the RPA correction 1 + r]Jj q Xo dominates the structure 
of the renormalizcd spin susceptibility Im^ and shows up 
as a strong commensurate resonance peak, although the 
bare spin susceptibility Im^o is incommensurate. When 
the wavevector q moves away from Q, the p-h transitions 
will move to the node direction and that connecting (0, 7r) 
and (7r,0), respectively. As a result, the single step-like 
rise splits into two as shown in Fig. 3(b). One is at the en- 
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ergy below it and the other above it, which correspond to 
the excitations E — > F and G — > H in Fig.l, respectively. 
The pole equation still satisfies around the high-energy 
peak (The line denoting —l/r)J~f q for q ^ Q will rise 
above the thin line in Fig. 3(a) due to \j q \ < \jq\), but 
the pole position falls well above the gap and the resul- 
tant p-h excitations are overdamped. On the other hand, 
though the peak caused by the low-energy jump does 
not satisfy the pole equation, because its jump height is 
reduced due to the weakening of the vH singularity at 
the q points E and F which are away from (0, — tt) and 
(71-, 0), it still affects the structure of the spin response 
via the RPA correction because of the small damping of 
the excitations. We can see from Fig. 3(b) that, when 
the deviation of the wavevector q from Q increases, the 
corresponding low-energy step-like rise decreases. From 
Ref. 0, we also know that the incommensurability of 
Imxo increases with the decrease of frequency. So, the 
peak dispersion determined both from the RPA correc- 
tion and the bare spin susceptibility Imxo has a down- 
ward curvature. But, which one mainly determines the 
position of the incommensurate peak in the renormalized 
spin susceptibility is interesting. This issue is related 
to the argument which attributes the incommensurate 
magnetic response and the commensurate resonance to 
belong to parts of the same collective mode In 
the previous investigation ||, the incommensurate mag- 
netic response is explained solely based on the structure 
of Imxo, because of the consideration x ~ Xo- However, 
we find that the peak position of Imx in q space, i.e. the 
incommensurability, is mainly determined by the RPA 
correction l + ?7</7gXo and in turn by the low-energy step- 
like rise. This can be seen clearly from Fig. 3(b), where 
the dashed-dotted line, the dotted line and the dashed 
line denote the imaginary parts of the bare spin suscep- 
tibility Imxo at q — (tt, 0.727t), (tt, 0.8tt) and (7r,0.97r), 
respectively. These wavevectors are the incommensu- 
rate peak positions in the renormalized spin suscepti- 
bility Imx at frequencies w — 0.15 J, 0.3 J and 0.45 J, 
respectively. Clearly, the low-energy step-like rises in 
Imxo shown in Fig. 3(b) and consequently the low-energy 
peaks in their real parts shown in Fig. 3(a) are at the 
same values as the above frequencies for the correspond- 
ing wavevectors. On the other hand, the incommensu- 
rabilities in Imxo for ui — 0.1 5 J, 0.3 J and 0.45 J are at 
q = (tt, 0.727r), (tt, 0.767r) and (tt, 0.887r). The same situa- 
tion is found for other doping concentrations. Therefore, 
both the resonance peak and the incommensurability are 
mainly determined by the step- like rise of Imxo- 

Below the lowest frequency for each doping density 
shown in Fig. 2, we find that a crossover of the incom- 
mensurate peaks from the (tt, tt + Sir) direction to the di- 
agonal direction occurs. It is because the excitations will 
be limited to be along the diagonal direction(the node- 
to-node excitation), such as E — > F shown in Fig.l, when 
the excitation energy is below the crossover energy due 



to the energy conservation law and the <i-wave symme- 
try of the SC gap. It is a specular feature for the model 
based on the FS topology and may also act to distinguish 
between the stripe and the FS nesting models. Unfortu- 
nately, the crossover energy for the underdoped system 
is too low to be available for the experimental observa- 
tion due to the dim scattering intensity |T^ |. However, 
we can see from Fig. 2 that it increases with doping den- 
sity and is about 0.135 J, 0.175 J and 0.2 J for the doping 
x = 0.08,0.14 and 0.20, respectively. So, we expect that 
it may be in the experimental observable range for the 
overdoped cuprates. For example, the crossover energy 
for x = 0.20 is 0.2J « 26meV which is experimentally 
accessible according to the recent report jl2| . Above the 
resonance frequency, the incommensurate peak occurs 
again. In particular, wc find that only in the underdoped 
regime (x < 0.13) are the incommensurate peaks along 
the (tt, tt + 8tt) direction. They are along the diagonal 
direction in the optimally doped and overdoped regimes 
(x > 0.14). It is in fact reflects the doping dependence of 
the intensity ratio of the incommensurate peak along the 
(tt, tt + Stt) direction to that along the diagonal direction, 
which decreases with the increase of doping density when 
the frequency is above the resonance energy pl| . From 
Fig. 2, one can see that the incommensurability above the 
resonance frequency increases with frequency, so the dis- 
persion shows a upward curvature which is similar to that 
of a massive magnon in a disorder Neel state. 

We show the doping dependences of the incommensu- 
rability S in the (tt, tt + 6tt) direction in Fig. 4(a) for ener- 
gies w = 0.15 J, 0.175 J, 0.20 J, 0.25 J and 0.35 J. In the un- 
derdoped regime, the incommensurability increases with 
the increase of doping and exhibits a nearly linear doping 
dependence for all energies. In the high-energy range, 
the incommensurability decreases with frequency. For 
La2- a; Sr :E Cu04, it is found that the incommensurability 
is energy independent for Two < 15meV. The same be- 
havior is found here when the frequency is reduced to be 
below about oj = 0.20 J, as can be seen from Fig. 4(a). By 
a closer inspection of the doping dependence of the in- 
commensurability at low frequencies such as to = 0.175 J 
and 0.15 J, we find that the equality S — 2x holds in the 
doping range from x = 0.06 to x = 0.13. We note that 
the incommensurability <5 defined here is twice that used 
by experimentalists which is indexed in units of the re- 
ciproal lattice vectors. This equality has been observed 
at low frequencies in underdoped cuprates |^,^| and is 
believed to be explained only in the stripe model be- 
fore 0. Therefore, our result provides an alternative 
explanation based on the FS nesting. However, the in- 
commensurability increases continually in the whole dop- 
ing range, which is different from the experimental ob- 
servation that it saturates near the optimal doping and 
in the overdoped regime PJj^fl- We also show the dop- 
ing dependences of the resonance energy E r in Fig. 4(b). 
This result has been reported by one of us (Li) and his 
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coworkers in a smaller doping range before |18|. Simi- 
larly, a nearly linear doping dependence is found in the 
underdoped regime, then it saturates in the slightly over- 
doped regime and eventually decreases with the further 
increase of dopings. Due to the similar doping depen- 
dence in the underdoped regime, the ratio of E r to 6 
shows a linear behavior as can be seen in the inset of 
Fig. 4(b). These results are remarkably in agreement with 
experiments |[^,[l6j. To understand these behaviors, we 
show the frequency dependences of the bare spin suscep- 
tibility xo a t the incommensurate q-point (jr, 0.847r) for 
different dopings in Figs. 3(c) and 3(d), and at the com- 
mensurate q-point (ir, ir) in the inset of Fig.3(d). We find 
that the step-like rise at (ir, tt) and the low-energy step- 
like rise at (tt, 0.847r) increases monotonously with doping 
in the underdoped regime. But the former saturates in 
the optimally doped regime and then decreases with the 
further increase of dopings, and the latter increases with 
dopings in the whole doping range considered here. Ac- 
cording to the above discussion, it is this feature that 
leads to the special doping dependences for both cases. 
From this reasoning, we may also argue that the satura- 
tion of the incommensurability observed experimentally 
near and after the optimal doping should not be related 
only to the change of the FS topology. Whether new 
physics such as the stripe phase is required or some ad- 
ditional factors such as the change of the quasiparticle 
nature due to a quantum critical point near the optimal 
doping [ p3[ should be taken into consideration is a open 
question. 

In summary, we have examined the doping and fre- 
quency dependences of the incommensurability and the 
resonance energy based on the Fermi surface topology. 
The calculations are carried out by use of the slave-boson 
mean-field approach to the t—t' — J model and includ- 
ing the antiferromagnetic fluctuation correction in the 
random-phase approximation. We find that the equality 
of the incommensurability and the doping density exists 
at lower frequencies in the underdoped regime, which has 
not been explained alternatively other than the stripe 
model before. The downward dispersion is reproduced 
and its doping dependence is presented for further ex- 
perimental testing. We also find a good linear behavior 
between the incommensurability and the resonance en- 
ergy in the underdoped regime, which is consistent with 
experiments. Our results may suggest a common origin 
for the incommensurate spin response and the resonance 
peak based on the Fermi surface topology and the <i-wave 
superconductivity. 
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FIGURE CAPTIONS 



Fig.l Fermi surface for the dispersion ek (thick solid 
line) at hole doping x = 0.12. The thin solid and dashed 
lines with arrow denote the threshold particle-hole exci- 
tations for different wavevectors(see text). 

Fig. 2 Peak positions of the imaginary part of the renor- 
malized spin susceptibility in the q space as a function 
of frequency for doping concentrations x = 0.08,0.14 
and 0.20. In the downward dispersion range, the in- 
commensurate peaks are along the (tt 7 it + Sir) direc- 
tion, so the horizontal axis is in unit of (l,q)n. Below 
ui = 0.135 J, 0.175 J and 0.2 J at dopings x = 0.08,0.14 
and 0.20, respectively, and for frequencies above the res- 
onance energy at x = 0.14 and 0.20, they are along the 
diagonal direction. In the latter case, the horizontal axis 
is in unit of (q, q)ir. 

Fig. 3 Frequency dependence of the bare spin suscep- 
tibility xo- Figures (a) and (b) are the real and imag- 
inary parts of the bare spin susceptibility at doping 
x=0.14 for different q. The dashed-dotted line, the dot- 
ted line and the dashed line denote spin susceptibilities 
at q = (ir, 0.727r), (it, 0.8ir) and (it, 0.97r), respectively. 
The thin line in figure (a) denotes —l/2i]J, so its cross 
with the real part corresponds to the pole of the renor- 
malized spin susceptibility. Figures (c) and (d) are the 
real and imaginary parts of the bare spin susceptibil- 
ity at q = (it, 0.847r) for different dopings. The solid 
line, dashed line and dotted line denote spin suscepti- 
bility at x — 0.08,0.14 and 0.2, respectively. The inset 
of figure (d) shows the frequency dependence of Im^o a t 
q = (it, 7t) for the doping x = 0.08(solid line),0.14(dashcd 
linc),0.20(dotted line) and 0.24 (dashed-dotted line). 

Fig. 4 (a) Doping dependences of the incommensurabil- 
ity S in the (it, tt + Sir) direction for different frequencies, 
(b) Doping dependence of the resonance energy E r . The 
inset of figure (b) shows the ratio between the incommen- 
surability and the resonance energy. The solid line in the 
inset is a guide to the eye. Note that the incommensura- 
bility 5 defined here is twice that used by experimentists 
which is indexed in units of the reciproal lattice vectors. 
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Fig. 1 Li et al. 
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Fig. 2 Lietal. 




Fig. 3 Li et al. 
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Fig. 4 Li et al. 



